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The critical behavior of semi-infinite systems in fixed di- 
mensions d < 4 is investigated theoretically. The appro- 
priate extension of Parisi's massive field theory approach is 
presented. Two-loop calculations and subsequent Pade-Borel 
analyses of surface critical exponents of the special and ordi- 
nary phase transitions yield estimates in reasonable agree- 
ment with recent Monte Carlo results. This includes the 
crossover exponent $(d = 3), for which we obtain the val- 
ues $(n = 1) ~ 0.54 and $(n = 0) ~ 0.52, considerably lower 
than the previous e-expansion estimates. 
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Methods of field theory have become an essential ingre- 
dient of the modern theory of critical phenomena, pro- 
viding both a conceptually appealing theoretical frame- 
work as well as powerful tools for quantitative analyses. 
Many important field-theoretical techniques, developed 
and tested originally in the study of bulk critical behavior, 
subsequently have been extended to systems with bound- 
aries and utilized with considerable success in the anal- 
ysis of surface critical phenomena at bulk critical points 

0. 

One important method that, to our knowledge, has not 
yet been employed in the study of surface critical phe- 
nomena is the field-theoretic renormalization group (RG) 
approach in a fixed space dimension d below the upper 
critical dimension d* ||. The merits of this approach 
are well known: pushed to sufficiently high orders of per- 
turbation theory and combined with sophisticated series 
resummation techniques, it has produced |||5| values of 
bulk critical exponents belonging to the most accurate 
ones obtained so far |@-@]. In addition, it has proved 
very useful in quantitative investigations of preasymp- 
totic critical behavior |P-[r^|. 

In the present Letter we show how this approach can 
be generalized to systems with surfaces. Aside from the 
fundamental importance of such an extension, we have 
been motivated by the discrepancies existing between the 
results of recent Monte Carlo analyses for the crossover 
surface exponent <& in three-dimensional systems 13 1(| 
and previous numerical estimates of $ obtained by ex- 
trapolation of its second-order e = 4 — d expansion to 
e = 1 |. In particular, the most recent Monte Carlo 
values for $ of both the Ising Jl5| and polymer |l6| uni- 
versality classes are substantially lower than the original 



e-expansion estimates, and also significantly lower than 
previous results of computer simulations |17 T^JT^]. The 
annoyingly large discrepancy with the e-expansion results 
calls for clarification and improvement of field-theoretic 
analyses. With this in mind, we have extended the mas- 
sive field theory approach || to the study of surface crit- 
ical behavior. Using this approach, we have performed 
two-loop calculations for the semi-infinite <^ =3 model. 
Subsequent Pade-Borel analyses of the resulting series 
expansions yield the estimates for surface critical expo- 
nents shown in Table 1. 

The Hamiltonian of the model is 



H= ( d d z 
Jv 



l d - l x^ca <t> 2 , 
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where <p = (feW) is an n-component field on the half- 
space V = {x = (x||,z) | X|| 6 M d ^ x ,z > 0} bounded 
by dV, the z = plane. For the time being we assume 
the regularization of parallel-momentum integrations by 
a large-momentum cutoff A. 

In wishing to study critical systems in fixed dimen- 
sions d < d* by means of perturbative methods, one is 
faced with the familiar occurrence of infrared (ir) sin- 
gularities. Feynman integrals involving massless propa- 
gators become ill-defined. This is a problem typical of 
massless sMperrenormalizable field theories. Nevertheless 
such theories exist, owing to the appearance of a non- 
perturbative critical bare mass shift |p|, ^8|,^9[ . In the di- 
mensionally regularized theory, this shift corresponds to 
a critical value of tq of the form tq c — u^ e T(e) , where 
T(e) has simple poles at e = e& = 2/fc, k £ IN. 

In the case of our semi-infinite model (Q), an inde- 
pendent mass scale is provided by cq, the bare surface 
enhancement. Its physical significance is well-known [0: 
depending on whether Cq is smaller, equal, or larger than 
a certain critical value c sp = c sp (e, Mo,A), distinct types 
of surface transitions occur at the bulk critical point, 
called ordinary, special, and extraordinary transitions, 
respectively. Finite values of c sp exist, of course, only if 
d exceeds the lower critical dimension d*(n) = 3 — S n \ for 
the occurrence of surface ordering. 

The quantity c sp has a number of properties analogous 
to tq c . First, in the cut-off regularized theory, it diverges 
as A — > oo. Its leading uv singularity for d < 4 is of the 
form uq A d ~ 3 , similar as tq c has one ~ uq A d 
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A second analogy becomes apparent if we consider 
the uv behavior of correlation functions. In the case 
of the 4>d<4 bulk theory (with V = M d ), a mass shift 
T o = T 0c + Stq is known to be sufficient to absorb the uv 
singularities of the A-point vertex and correlation func- 
tions. In our semi-infinite case we must keep track in cor- 
relation functions of which points are taken on or off the 
surface. Let G^ N ' M '({xi}, {r.,}; r , Uq, cq, A) be the (reg- 
ularized) connected (N + M)-point correlation functions 
involving N fields </>(Xj) at distinct points x$, 1 < i < N, 
off the surface and M fields <j>{vj ,z = 0) = cj) s (r 3 ) at dis- 
tinct surface points with parallel coordinates r j , 1 < j < 
M. To absorb the uv singularities of these functions for 
d < 4 both a mass shift and a surface-enhancement shift 
c o = c sp + <5co is required; i.e., the limits A — > oo of the 
G<- N < M ) with fixed 5t q > and 6c > exist. 

A further analogy is the form of c sp in the dimension- 
ally regularized theory. On dimensional grounds one ex- 
pects it to read 



C(e), 



(2) 



where C(e) has poles at the above values The value 
of c sp is fixed by the requirement that the layer sus- 
ceptibility Xu = G( 0,2 )(p = 0) diverges at the special 
point (roc,c sp ). Here & ' 2 \p) is the parallel Fourier 
transform of G^ 0,2 '(r — r'). Hence c sp is the solution to 
Xii (to c , uq, Cq = c sp ) = . By explicit computation of 
the 0(uq) term in the expansion of Xii one may verify 
that it has indeed a pole at d = 3, which must be elim- 
inated by a surface-enhancement shift. Thus, just as in 
the bulk theory, one encounters the following situation: 
although, by making appropriate shifts of To and cq, the 
theory can be made uv-finite for d < 4, the critical values 
TOc and c sp cannot be determined by perturbation theory 
owing to their nonanalytic dependence on Uq. In order 
to obtain meaningful results in spatial dimensions d with 
d* < d < d* one should employ an appropriate extended 
version of the massive field theory RG approach, avoiding 
direct computations of both tq c and c sp . 

To explain our method, we start by recalling from Ref. 
jjj that the field 0(x) and the surface field 4> s (r) should be 
reparametrized by distinct renormalization factors. Ac- 
cordingly we introduce the renormalized correlation func- 
tions through 

GgM ( .; m, u, c) = Z - {N+M)/2 Z~ M/2 G^ . (3) 

Unlike the renormalized surface enhancement c and 
the renormalization factor Z\ (which are specific to 
our surface-bounded system), the renormalized mass m, 
renormalized coupling constant u, and Z^ are bulk quan- 
tities. We fix the latter ones by standard normaliza- 
tion conditions for the renormalized bulk vertex functions 
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In terms of their Fourier transforms V 
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In order to specify c and Z\, we require that 

1 



G( c ° n 2) (p;™,^c) 



p=0 
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p 2 =o 2m(m + c) 2 
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These normalization conditions are suggested by the fa- 
miliar zero-loop expression |Q [c + (p 2 + To) 1 / 2 ] -1 for 
G (0 ' 2) (p). Together with (g), the first one, (0), ensures 
that the special point is located at m = c = 0. Equation 
(0) fixes Z\ in much the same way as (||) determines Z$. 
Since we are also interested in the crossover exponent 
we should also consider correlation functions with in- 
sertions of the surface operator (j) 2 . Let G^ N,M;1,11 ^ de- 
note the connected correlation function with N fields 4> 
off the surface, M surface fields <j> s% I insertions of \§ 2 
at points off the surface, and 1\ insertions of \<\> 2 S - We 
write the renormalized operators as (</> 2 ) re n = 4 >2 and 
(^s)ren = Z c cf) 2 . The bulk factor Z^i we fix through 

(2,1)/ 



the usual normalization condition T b ' R ({0}) = 1 for the 
Fourier transform of the bulk two-point vertex function 
. To specify Z c we require that 



with an insertion of 



^ 2! °' 1) ({p};m,«,c)| {p=0} = (m + c)- 



(9) 



these conditions read 



Again, this condition is suggested by the zero-loop result 
for G (0 < 2) , as can be seen by comparison with y7|) , noting 
that the bare analog of the left-hand side of (^ may be 
written as (—d/dco) G (0 ' 2 )(0). 

The above normalization conditions determine r , u , 
Z$, Z^t, Co, Zi, and Z c as functions of u, c, and A. 
All Z-factors have a finite A — > oo limit in the d < 4 
case considered here. Let us set A = oo in the sequel. 
(In our calculations we actually took A = oo from the 
outset, employing dimensional regularization.) Then the 
bulk Z-factors Z$ and Z^i become functions of the single 
dimensionless variable u. However, with our choice of 
normalization conditions, the surface Z-factors Z\ and 
Z c depend on both u and the dimensionless ratio c/m. 

By varying m at fixed uq, cq, and A, the analogs of 
the Callan-Symanzik equations (CSE) can be derived in 
a straightforward fashion. Let us define the functions 
(3(u) = md m \ u, B c (u,c/m) = md m \ Q hxc, r)$(u) = 
md m \o lnZ^, and r)x(u,c/m) = md m \o InZi, where |o 
indicates that the derivatives are taken at fixed uq, Co, 
and A. Upon introducing the differential operator 
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P m = 77i d m + fid u + B c cd c ^ — 770 + — rji 



(10) 



the CSE of the Gi^l can be written as 



V m G(f n ' M) = -(2 - ^) m 2 / d d X GW-°) , (11) 



where the integration is over the position X of the in- 
serted (j) 2 operator. 

Just as in the bulk case, and as could be corroborated 
by means of a short-distance expansion, the inhomogene- 
ity on the right-hand side should be negligible in the crit- 
ical regime. The resulting homogeneous CSE differs from 
its standard bulk analog in that it involves functions /3 C 
and r]i of two variables, u and c/m. In a study of the 
crossover from the critical behavior characteristic of the 
special transition (for c/m <C 1) to that of the ordinary 
transition (for c/m 3> 1) it would be essential to carry 
along this dependence on c/m. Such a crossover analysis 
is beyond the scope of this Letter. Our main aim here is 
the calculation of the surface critical exponents of these 
two surface transitions. Therefore we focus directly on 
the respective asymptotic behavior. 

The easier case is the special transition, which we con- 
sider first. Its critical behavior can be investigated upon 
setting c = 0. The homogeneous CSE can then be inte- 
grated in a standard fashion. One finds that the surface 
correlation exponent can be written as 



sp sp/ *\ . 



(12) 



where r\ = 770(77*) is the standard bulk correlation ex- 
ponent, u* is the ir-stable zero of f3(u), and tjI p (u) = 
771 (u, 0). A similar analysis applied to the CSE for 
correlation functions with insertions of (j) 2 s shows that 
the anomalous dimension of this operator is given by 
the fixed-point value of the function r) c (u,c/m) = 
md m \o lnZ c . Hence the surface crossover exponent can 
be expressed as 



$ = v[l + r)l p (u*)] 



(13) 



in terms of ti^P(u) = rj c (u, 0) and the bulk exponent v. 

We have computed the perturbation expansions of the 
functions "ff^(u) = r\^(u) + 77^ (u) and ill P (u) to two- 
loop order for d — 3. Let us make the usual change of 
normalization of u such that (3(u) = —u + u 2 + 0(u 3 ). 
Then our results become 



n + 2 12 (n + 2) 

"11 {u) = ~^T8) U 



(n + i 



2A 
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and 



Z±* (1*4-1 



24 (n + 2) 



(n + 8) 2 
n + 2 



A-B 



In 2 2- T ln2+1? 



12 

1 

32 
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The numbers A and B are associated with contributions 
of two-loop diagram with three equivalent internal lines. 
We have determined them by numerical integration, ob- 
taining A = 0.202428 and B = 0.678061. 

From these results analogous series expansions can be 
derived for other surface exponents by means of scaling 
laws p]. The coefficients of most of the resulting series 
alternate in sign and decrease in absolute value |Q . Ex- 
ceptions are some series involving 77^, whose behavior is 
rather bad. We have performed Pade analyses and Pade- 
Borel resummations of the direct and inverse series for 
the indices 77||, Ai, ?7_l, Pi, 711, 71, Si, Su, eta, an, and $. 
In these computations we have used the fixed-point val- 
ues u*(n = 0) = 1.632 and u*(n = 1) = 1.597 from the 
Pade-Borel calculations of Refs. pl[ and ps| . 

In the group of exponents 7711 , . . . , Su related to 77^ (u*) 
the most reliable estimate is obtained for Ai, which ap- 
pears to exhibit the best convergence properties. Pade- 
Borel resummations of this series yield the values A 1(77 = 
0) = 0.921 and A 1 (n = 1) = 0.997. The estimates in Ta- 
ble 1 for the other exponents of this group have been 
derived from these values of Ai via scaling relations, us- 
ing the d = 3 values of v and 77 given in Ref . || . 

In the group of exponents cti, an, $ related to r)^?(u*) 
the series for ai yields estimates with the least scatter- 
ing. The values of an and $ listed in Table 1 have been 
obtained from our results for a\ and the accepted bulk 
values of v and 77. 

In order to study the ordinary transition, we must con- 
sider the limit 777 —> at fixed c > 0; i.e., c/777 — > 00. A 
well-known complication is that one cannot simply set 
c = 00 in the functions G^ N ' M ^ with M > because 
a Dirichlet boundary condition holds for c = 00 [0. 
Yet, the cumbersome calculation of c-dependent quan- 
tities can be avoided by considering correlation functions 
involving instead of <j) s the normal derivative d n (f> at the 
surface. Let be the analog of G^ N,M ^ resulting by 

replacement of all the <f> s by d n <f>, with Co = 00. Following 
the strategy described in Ref. HI], we have performed an 
independent RG analysis of these functions. The renor- 
malized operator {d n 4>)n = (Z^Zi^) -1 / 2 d n 4> involves a 
renormalization factor Zi l00 (u), which we fix through the 
condition (d/dp 2 )G ( °'%(p,m,u)\ p2=0 = -I /2m |o||§. 

The associated RG function r}i j00 (u) = m d m \o In Zi iOC 
gives us the function 

77f' d (7i) = 2 + 771,00 (u) +v4 u ) (16) 

needed to determine the surface exponent 77|| rd = 



3 



rf^ ld (u*). The explicit two- loop expression of this func- 
tion for d — 3 reads 



24 (n + 2) 



C 



14 



(17) 



with C = -0.105063. 

Again, this result can be combined with scaling laws to 
obtain the corresponding expansions of the other surface 
critical indices of the ordinary transition. In most cases 
the coefficients of these series do not alternate in sign, 
with similar behavior of the inverse ones. Therefore, they 
are not adapted to Pade-Borel resummation. The numer- 
ical values of the ordinary-transition exponents arising 
from [1/1] Pade approximants for each individual series 
are listed in Table 1. The exception is 7°£ d for which this 
approximant does not exist. We evaluated this index us- 
ing the values of Tyj? and v. For n — 2 and n — 3, the 
fixed-point values u* = 1.558 and u* = 1.521 [^2j were 
used, respectively. 

Our numerical values of the surface critical exponents 
gathered in Table 1 generally are in reasonable agreement 
both with previous estimates based on the e expansion 0] 
as well as with those obtained by other means (see Refs. 
and fl3| - |l7| , p4f for comparisons). Note, however, that 
our estimates for $ are definitively lower than the values 
$(n = 1) ~ 0.68 and = 0) ~ 0.67 quoted in Ref. §, 
which were obtained by setting e = 1 in the e expansion 
of $ to order e 2 . Recent Monte Carlo simulations have 
yielded the significantly lower estimates $(1) = 0.461 ± 
0.015 pi, $(0) = 0.530 ±0.007 @, and $(0) = 0.496 ± 
0.005 |l§. Our results $(1) ~ 0.54 and $(0) ~ 0.52 
are fairly close to these values. This indicates that the 
crossover exponents $ for d = 3 and n = 0, 1 are indeed 
smaller than previously thought. Let us also note that 
when the e expansion to order e 2 is extrapolated to d = 3 
by means of more elaborate techniques (e.g, Pade-Borel 
resummation), estimates for $(0) and $(1) as low as ours 
can be obtained. 

A more detailed account of our work will be given else- 
where [ po[ . 
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TABLE I. Estimated d = 3 values of the surface critical 
exponents of the special and ordinary transition for various 
values of n. 





SPECIAL 




ORDINARY 






n = 


n = 1 


n = 


n = 1 


n = 2 


n = 3 


V\\ 


-0.133 


-0.165 


1.660 


1.528 


1.422 


1.338 


Ai 


0.921 


0.997 


0.382 


0.450 


0.514 


0.574 


r?± 


-0.053 


-0.067 


0.859 


0.801 


0.753 


0.714 


ft 


0.255 


0.263 


0.817 


0.845 


0.868 


0.889 


7n 


0.666 


0.734 


-0.388 


-0.333 


-0.282 


-0.238 


7i 


1.207 


1.302 


0.664 


0.742 


0.815 


0.882 


Si 


5.734 


5.951 


1.854 


1.941 


2.019 


2.088 


5n 


3.612 


3.791 


0.504 


0.582 


0.651 


0.711 


ai 


0.342 


0.279 










an 


-0.140 


-0.182 










$ 


0.518 


0.539 
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